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ABSTRACT 

There is a growing evidence that extended radio halos are most likely generated 
by electrons reaccelerated via some kind of turbulence generated in the cluster volume 
during major mergers. It is well known that Alfven waves channel most of their en- 
ergy flux in the acceleration of relativistic particles. Much work has been done recently 
to study this phenomenon and its consequences for the explanation of the observed 
non-thermal phenomena in clusters of galaxies. We investigate here the problem of 
particle-wave interactions in the most general situation in which relativistic electrons, 
thermal protons and relativistic protons exist within the cluster volume. The inter- 
action of all these components with the waves, as well as the turbulent cascading 
and damping processes of Alfven waves, are treated in a fully time-dependent way. 
This allows us to calculate the spectra of electrons, protons and waves at any fixed 
time. The Lighthill mechanism is invoked to couple the fluid turbulence, supposedly 
injected during cluster mergers, to MHD turbulence. We find that present observations 
of non-thermal radiation from clusters of galaxies are well described within this ap- 
proach, provided the fraction of relativistic hadrons in the intracluster medium (ICM) 
is smaller than 5 — 10%. 

Key words: acceleration of particles - radiation mechanisms: non-thermal - galaxies: 
clusters: general - radio continuum: general - X-rays: general 



1 INTRODUCTION 

There is now firm evidence that the ICM is a mixture of hot 
gas, magnetic fields and relativistic particles. While the hot 
gas results in thermal bremsstrahlung X-ray emission, rela- 
tivistic electrons generate non-thermal radio and hard X-ray 
radiation. The amount of the energy budget of the intraclus- 
ter medium in the form of high energy hadrons can be large, 
due to the phenomenon of confinement of cosmic rays over 
cosmological time scales (Volk et al. 1996; Berezinsky, Blasi 
& Ptuskin 1997). Nevertheless, the gamma radiation that 
would allow us to infer the fraction of relativistic hadrons in 
clusters has not been detected as yet (Reimer et al., 2003). 

The most important evidence for relativistic electrons in 
clusters of galaxies comes from the diffuse synchrotron radio 
emission observed in about 35% of the clusters selected with 
X-ray luminosity > 10 45 erg s" 1 (e.g., Feretti, 2003). The 
diffuse emission comes in two flavors, referred to as radio 
halos (and/or radio mini-halos) when the emission appears 



concentrated at the center of the cluster, and radio relics 
when the emission comes from the peripherical regions of 
the cluster. 

The difficulty in explaining the extended radio halos 
arises from the combination of their ~Mpc size, and the rela- 
tively short radiative lifetime of the radio emitting electrons. 
Indeed, the diffusion time necessary for the radio electrons to 
cover such distances is orders of magnitude larger than their 
radiative lifetime. As proposed first by Jaffe (1977), a solu- 
tion to this puzzle would be provided by continuous in situ 
reacceleration of the relativistic electrons on their way out. 
This possibility was studied more quantitatively by Schlick- 
eiser et al. (1987) who successfully reproduced the integrated 
radio spectrum of the radio halo in the Coma cluster. In the 
framework of the in situ reacceleration model, Harris et al. 
(1980) first suggested that cluster mergers might provide the 
energetics necessary to reaccelerate the relativistic particles. 

An alternative to the reacceleration scenario was put 
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forward by Dennison (1980), who suggested that relativis- 
ts electrons may be produced in situ by inelastic proton- 
proton collisions through production and decay of charged 
pions. This model is known in the literature as the secondary 
electron model. 

Diffuse radio emission is not the only evidence of non- 
thermal activity in the ICM. Additional evidence, although 
limited to a few cases, comes from the detection of extreme 
ultra-violet (EUV) excess emission (e.g., Bowyer et al. 1996; 
Lieu et al., 1996; Berghofer et al., 2000; Bonamente et al., 
2001), and of hard X-ray (HXR) excess emission in the 
case of the Coma cluster, A2256 and possibly A754 (Fusco- 
Femiano et al., 1999, 2000, 2003, 2004; Rephaeli et al., 1999; 
Rephaeli & Gruber, 2002, 2003)*. While, with the exception 
of the Coma and Virgo clusters, the detection of EUV emis- 
sion is still controversial (e.g., Berghofer et al., 2000), the 
detection of HXRs appears robust as it is claimed indepen- 
dently by different groups and with different X-ray observa- 
tories (BeppoSAX, RXTE). If these excesses are indeed of 
non-thermal origin, they may be explained in terms of IC 
scattering of relativistic electrons off the photons of the cos- 
mic microwave background (CMB) (Fusco-Femiano et al., 
1999, 2000; Rephaeli et al., 1999; Volk & Atoyan 1999; Blasi 
2001; Brunetti et al.2001a; Fujita & Sarazin 2001; Petrosian 
2001; Kuo et al., 2003). Alternatively HXRs might also re- 
sult from bremsstrahlung emission of electrons in the ther- 
mal gas whose spectrum is not a perfect Maxwellian, but 
rather has a supra-thermal tail (e.g., Ensslin, Lieu, Bier- 
mann 1999; Blasi 2000; Dogiel 2000; Sarazin & Kempner 
2000). Blasi (2000) calculated the spectrum of the elec- 
trons as modified by the resonant interaction with MHD 
waves by solving a time-dependent Fokker-Planck equation 
and reached the conclusion that the intracluster medium is 
heated and at the same time develops a supra-thermal tail 
due to the resonant particle-wave interaction. An energy in- 
jection comparable with the whole luminosity of a merger, 
lasting for about a billion years was required in this calcula- 
tion (Blasi 2000) in order to explain observations. A shorter 
allowed duration (about 10 8 years) of the event was esti- 
mated by Petrosian (2001). 

Both the IC model and the bremsstrahlung interpre- 
tation have problems: the first one would require cluster 
magnetic field strengths smaller than those inferred from 
measurements of the Faraday rotations (e.g., Carilli & Tay- 
lor 2002). The second one would require a large amount of 
energy to maintain a substantial fraction of the thermal elec- 
trons far from the thermal equilibrium for more than a few 
10 8 yrs (Blasi, 2000; Petrosian, 2001). 

The origin of radio halos is still an open problem. In 
principle, it is known that the fine radio properties of the 
radio halos may be naturally accounted for in models in 
which electrons are directly accelerated (primary electrons) 
whereas they are difficult to reproduce in models in which 
the electrons are secondary products of hadronic interac- 
tions (Brunetti, 2002). On the other hand, the details of the 
physics used in the primary electron models are relatively 
uncertain. In principle, merger shocks can accelerate rela- 
tivistic electrons to produce large scale synchrotron radio 



* The evidence for HXR excess in the case of the Coma cluster 
has been recently disputed by Rossetti & Molendi (2004) 



emission (e.g., Roettiger et al., 1999; Sarazin 1999; Takizawa 
& Naito, 2000), however, the radiative life-time of the emit- 
ting electrons diffusing away from these shocks is so short 
that they would just be able to produce relics and not Mpc 
scale radio halos (e.g., Miniati et al, 2001). In addition, a 
number of papers (Gabici & Blasi 2003; Berrington & Der- 
mer 2003) have recently pointed out that the Mach number 
of the typical shocks produced during major merger events 
is too low to generate non-thermal radiation with the ob- 
served fluxes and spectra. 

Re-acceleration of a population of relic electrons by tur- 
bulence powered by major mergers is suitable to explain 
the very large scale of the observed radio emission and is 
also a promising possibility to account for the fine radio 
structure of the diffuse emission (Brunetti et al., 2001a,b). 
There are a number of possibilities to channel the energy of 
the turbulence in the acceleration of fast particles, namely 
via Magneto-Sonic (MS) waves, via magnetic Landau damp- 
ing (e.g., Kulsrud & Ferrari 1971), via Lower Hybrid (LH) 
waves (e.g., Eilek & Weatherall 1999) or via Alfven waves. 
Since Alfven waves are likely to be able to transfer most 
of their energy into relativistic particles, they have received 
much attention in the last few years. In this framework for 
instance Ohno, Takizawa and Shibata (2002) developed a 
time-independent model for the acceleration of the relativis- 
tic electrons expected in radio halos through magnetic tur- 
bulence. The authors studied the acceleration of continu- 
ously injected relativistic electrons by Alfven waves with a 
power law spectrum and applied this model to the case of 
the radio halo in the Coma cluster. More recently, Fujita, 
Takizawa and Sarazin (2003) studied the effect of Alfvenic 
acceleration of relativistic electrons in clusters of galaxies. 
These authors invoked the Lighthill theory to establish a 
connection between the large scale fluid turbulence and the 
radiated MHD waves. The electron and MHD-wave spectra 
adopted by Fujita et al.(2003) are obtained via a self-similar 
approach by requiring that the spectra are described by two 
power laws. 

These approaches have two intrinsic limitations: the 
first one is in the assumption, mentioned above, that all 
spectra are time-independent and that the turbulence spec- 
trum is a power law. The second is that they neglect, as 
all other previous approaches did, the effect of relativistic 
hadrons in the ICM: it is well known that the interaction 
of the Alfven waves with relativistic particles is, in gen- 
eral, more effective for protons than for electrons (e.g., Eilek 
1979). It is also well known that the presence of a significant 
energy budget in the form of relativistic particles can sig- 
nificantly affect the spectrum of the Alfven waves through 
damping. In fact, this damping occurs even on the thermal 
protons in the ICM, another effect which was never included 
in previous calculations. 

The calculations presented here provide a self-consistent 
time-dependent treatment of the non-linear coupling of 
Alfven waves, relativistic electrons, thermal and relativis- 
tic protons. The results previously appeared in the litera- 
ture can be obtained as special cases of our very general 
approach, which is in principle applicable to scenarios other 
than clusters of galaxies. 

The paper is organized as follows: in Sect. 2 we discuss 
the energy losses of relativistic particles and the presence of 
these particles in galaxy clusters. In Sect. 3 we discuss the 
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physics of Alfven waves, their generation in galaxy clusters 
and their interaction with particles. In Sect. 4 we discuss the 
physics of the coupling between Alfven waves, leptons and 
hadrons and derive the time evolution of waves and particles 
as a function of the physical conditions typical of the ICM. In 
Sect. 5 we apply our general formalism to calculate the fluxes 
of radio radiation and hard X-ray tails in galaxy clusters. 

Throughout the paper we adopt Ho = 50 km s _1 
Mpc -1 ; if not specified all the quantities are given in c.g.s. 
units. 



2 BASIC EQUATIONS AND ASSUMPTIONS 

2.1 Energy losses of relativistic particles in the 
ICM 

In this Section we give a short summary of the main en- 
ergy loss channels that may be important for electrons and 
protons in the ICM. 



2.1.1 Electrons 

Relativistic electrons with momentum p = m e crf in the ICM 
lose energy through ionization losses and Coulomb collisions 
(Sarazin 1999): 



\dtji 



-3.3 x 10"^n th 



1 + 



ln(7/n t h) 
75 



(1) 



where nth is the number density of the thermal plasma. Rel- 
ativistic electrons lose energy via synchrotron emission and 
inverse Compton scattering (ICS): 
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where B^a is the magnetic field strength in fj,G and 6 is 
the pitch angle of the emitting electrons; in case of efficient 
isotropization of the electron momenta it is possible to aver- 
age over all possible pitch angles, so that < sin 2 6 >= 2/3. 
It is well known that in the typical conditions of the ICM 
radiation losses are the most important for electrons with 
Lorentz factor 7 ^> 100 while Coulomb losses dominate at 
lower energies (Sarazin 1999,2001; Brunetti 2002). The life- 
time of relativistic electrons, defined as r ~ 7/7, can be 
easily estimated from Eqs.(l-2) as: 



1.0 >J (300) 
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2.1.2 Protons 



The main channel of energy losses for relativistic protons 
is represented by inelastic proton-proton collisions. The 
timescale associated with this process is : 



1 



10 1 



n t h 



(4) 

nthO-ppC VlO S J 

Inelastic pp scattering is weak enough to allow for the ac- 
cumulation of protons over cosmological times (Berezinsky, 
Blasi & Ptuskin 1997). The rare interactions with the ICM 



may generate an appreciable flux of gamma rays and neutri- 
nos, in addition to a population of secondary electrons (Blasi 
& Colafrancesco 1999). The process of pion production in pp 
scattering is a threshold reaction that requires protons with 
kinetic energy larger than ~ 300 MeV. 

Protons which are more energetic than the ther- 
mal electrons, namely protons with velocity j3 > f3 c — 
(3/2m e /m p ) 1 ^ 2 (3 e (/3 e here is the velocity of the thermal elec- 
trons, j3 e ~ 0.18(T/10 8 7f) 1/2 ) lose energy due to Coulomb 



interactions. If we define x r , 
(e.g., Schlickeiser, 2002): 
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with the following asymptotic behaviour: 



dp ( n th , 

d^viF^M p 
{ p 



Const, for p >> mc 

' for mcx m < p << mc (6) 



for mc(3 c < p < mcx Ti 



The timescale associated with Coulomb collisions (in the 
case mcxm. < p <S mc) can be therefore written as: 



t c ~ 2.5 x t„ 



\m v cj 



(7) 



For trans-relativistic and sub-relativistic protons this chan- 
nel can easily become the main channel of energy losses in 
the ICM. 



2.2 Origin and spectrum of the relic relativistic 
particles 

In this section we briefly discuss the mechanisms respon- 
sible for the injection of cosmic rays in galaxy clusters, 
more extended discussions have been previously presented 
by Berezinsky, Blasi & Ptuskin (1997), Kronberg (2002), 
Biermann et al. (2002) and Jones et al. (2002). 

Collisionless shocks are generally recognized as effi- 
cient particle accelerators through the so-called "diffusive 
shock acceleration" (DSA) process (Drury, 1983; Blandford 
& Eichler 1987). This mechanism has been invoked several 
times as the main acceleration process in clusters of galax- 
ies that have been involved in a merger event (Takizawa & 
Naito, 2000; Blasi 2001; Miniati et al., 2001; Fujita & Sarazin 
2001). For the case of standard newtonian shocks, relevant 
for clusters of galaxies, the spectrum of the accelerated par- 
ticles can be shown to be a power-law with a slope that 
is independent of the details of the diffusion in the shock 
vicinity, and depends only on the compression factor r at 
the shock, s = (r + 2)/(r - 1). 

The injection of CR ions at shocks is generally com- 
puted in the test particle limit while the injection efficiency 
is sometimes just assumed as a free parameter, while in other 
cases it is estimated according to the so-called thermal leak- 
age model (e.g., Kang & Jones, 1995). 

There is still some debate on the typical Mach number 
of the shocks developed in the ICM during cluster mergers. 
Some results from numerical simulations suggest the pres- 
ence of a large fraction of high Mach number shocks in clus- 
ter mergers (Miniati et al., 2000, 2001). Semi-analytical cal- 
culations (Gabici & Blasi 2003; Berrington & Dermer 2003) 
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Figure 1. Electron spectrum at z = injected as a single burst 
at Zi =0.01, 0.1, 0.3, 0.5 (from right to left) adopting injection 
spectrum Q(p) oc p~ 2,5 and maximum Lorentz factor ^/max — 
10 4 . The calculations arc carried out for n t ^ 



10" 



1/iG. 



have pointed out that the Mach numbers of the shocks re- 
lated to major mergers are expected to be of order unity. Re- 
cent numerical simulations (Ryu et al., 2003) seem to find 
more weak shocks than in Miniati et al.(2000, 2001). The 
comparison however with analytical calculations appears dif- 
ficult because of a different classification of the shocks in the 
two approaches. Moreover the numerical simulations of Ryu 
et al. (2003) also find a class of high Mach number shocks 
that are related to gravitationally unbound structures, not 
included in semi-analytical calculations. On the other hand, 
the weakness of the merger-related shocks seems to be also 
suggested by the few observations in which the Mach number 
of the shock can be measured (e.g., Markevitch et al., 2003). 
If shocks related to major mergers are indeed weak, the spec- 
tra of the accelerated particles are typically too steep to be 
relevant for nonthermal phenomena in clusters of galaxies. 

A relevant contribution to the injection of cosmic rays 
in clusters of galaxies may come from Active Galactic Nu- 
clei (AGN). AGNs indeed inject in the ICM a considerable 
amount of energy in relativistic particles and also in mag- 
netic fields, likely extracted from the accretion power of their 
central black hole (Ensslin et al., 1997). It should be stressed 
that the presence of relativistic plasma in AGNs (radio lobes 
and jets) is directly observed because of the synchrotron and 
inverse Compton emission of accelerated electrons. 

Finally, powerful Galactic Winds (GW) can inject rel- 
ativistic particles and magnetic fields in the ICM (Volk & 
Atoyan 1999). Although the present day level of starburst 
activity is low, it is expected that these winds were more 
powerful during starburst activity in early galaxies. Some 
evidence that powerful GW were more frequent in the past 
comes from the observed iron abundance in galaxy clusters 
(Volk et al. 1996). 



Since most of the scenarios discussed above imply power 
law spectra of relativistic particles at the injection sites, in 
the following we will restrict our calculation to this case. It 
is worth recalling that transport effects and energy losses 
modify the shape of these spectra, that are not expected to 
be power laws at later times. 

2.2.1 Electron Spectrum 

In the conditions typical of the ICM, ultra-relativistic elec- 
trons rapidly cool down through ICS and synchrotron emis- 
sion, and accumulate in the region of Lorentz factors 7 ~ 
100 — 500 where they may stay for a few billion years before 
cooling further down in energy through Coulomb scattering 
and eventually thermalize. 

The kinetic equation that describes these losses and the 
continuous injection of electrons is (Kardashev, 1962): 

dN e (p,t) _ d (.dp 



dt 



dp \ dt ra . 



+ -^)N c (p,t)\+Q c (p) (8) 
d dti J 

where dp/dt v and <ip/di rad are given by Eqs.(l) and (2) re- 
spectively and Q c (p) represents the injection term. The time 
evolution of relativistic electrons in the ICM during cosmo- 
logical times has been investigated in some detail by Sarazin 
(1999) by making use of the numerical solutions of Eq. (8). 
In Fig. 1 we plot the spectrum of relativistic electrons mea- 
sured at the present time and injected in a single event at 
z = Zi in the cluster volume. The calculations are shown for 
different values of Zi\ the spectra flatten with increasing Zi 
(typical conditions of the ICM are adopted, as described in 
the figure caption). Most electrons injected at z% > 0.2 get 
thermalized. 

In the case of continuous time-independent injection, it 
is well known that an equilibrium spectrum of the electrons 
is achieved, that can be written as follows: 



N e (p,z) oc i i{p>>peci{z) 



(9) 



where p cq is close to the momentum at which the Coulomb 
losses dominate over the radiative losses, namely: 

1/2 



Peci(z) ~ 300 x m c c 



IO- 3 



(1+*)- 



(10) 



More realistic injection histories of relativistic electrons 
can be easily implemented in this kind of calculation. 

2.2.2 Protons 

Both timescales of energy losses and diffusion out of the 
cluster volume are larger than the Hubble time for most 
of the cosmic ray protons in the ICM (Volk et al. 1996; 
Berezinsky, Blasi & Ptuskin 1997), although the energy at 
which confinement becomes inefficient is rather sensitive on 
the adopted diffusion model. 

As stressed above, mildly and sub- relativistic protons 
may be significantly affected by Coulomb energy losses, 
which in turn change the particle spectrum with respect to 
the injection spectrum. 

Below the threshold for pion production in the pp colli- 
sions (or simply by neglecting the effect due to pp collisions) , 
the time evolution of the proton spectrum is described by 
the following equation: 
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Log (p/m p c) 

Figure 2. Present-epoch spectrum of the cosmic ray protons con- 
tinuously injected (with p > Q.lm p c) in the ICM starting from 
z;. The spectra are plotted for s =2.4 (dotted lines) and 3.0 (solid 
lines) and for z; =0.1, 0.5, 1 (from bottom to top). Calculations 
are carried out assuming = 10~ 3 cm~ 3 . 
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Figure 4. Energy density in the form of cosmic ray protons in 
the ICM (arbitrary units) if protons are injected at a constant 
rate (number of particles per unit time) at p > 0.1 m p c with a 
power law injection spectrum with slope s (on the x-axis). The 
time at which injection starts is taken as Zj = 1. The density of 
the thermal gas is assumed to be n t ^ = 10 _3 cm~ 3 . 
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Figure 3. The ratio between the energy density of cosmic ray 
protons with momentum p > 0.1 m p c at the present time (z = 0) 
and the energy density in cosmic rays injected from the redshift 
Zi = z (x-axis) to the present. Calculations are carried out for 
s = 3.3 (solid line) and s = 4.0 (dashed line) and for n th = 
10 _3 cm- 3 . 



dN p (p,t) 
dt 



(ii) 



where dp/dt i is given by Eq. (5). Adopting the non- 
relativistic asymptotic behaviour of Eq. 6, the asymptotic 
solution (for time independent injection) is as follows: 



N p (p,T) = 



nth C 



[3 > x m 

4^{l-exp(^F)} 

where C is the constant in Eq. (5) and T is the time elapsed 
from the beginning of the injection. Eq. (12) presents the 
following behaviours : 



{p 3 3 if p « p»(T) & p > m p cx n 
p~ s , if p « p*(T) & p < nipCXv 
p- s , ifp»p*(T); 

where p, (T) is given by t : 



p,(T) ~ 0.7 x 



( T 



\ lOGyr 10" 



1/3 



(13) 



(14) 



In Fig. 2 we plot the spectrum of the protons measured 
at the present time as obtained solving Eq. (11) numerically 
under the assumption of a continuous time-independent in- 
jection of protons (starting from different Zi, see caption) 



t Note that in the ultra-relativistic limit Coulomb 
ligiblc and one has N p (p,t) tx p~ s 



losses are neg- 
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in the cluster volume. The proton spectra, largely modified 
in the trans-relativistic regime by Coulomb losses, has the 
asymptotic behaviour described by Eq. (13). 

Since the energy flux of Alfven waves, as we show below, 
is efficiently damped on protons, it is important to estimate 
the amount of energy stored in the form of protons in the 
ICM. 

In Fig. 3 we plot the ratio between the energy injected 
in cosmic ray protons (assuming time independent injection 
starting at redshift zi) and the energy stored in the form of 
supra-thermal protons (here p>0.1x m v c) at the present 
time for different values of the injected spectral index s (see 
caption). As expected, the ratio between the energy stored 
in protons at the present time and the injected energy is 
smaller for steeper injection spectra, due to the effect of 
Coulomb losses. 

Finally, by assuming that a fixed fraction of the ther- 
mal protons in the cluster volume is accelerated to higher 
energies at a constant rate (starting from Zi = 1.0), in Fig. 4 
we plot the energy stored in the cluster at the present time 
as a function of the slope of the injected spectra s. 



3 ALFVENIC ACCELERATION OF 
RELATIVISTIC PARTICLES 

Alfven waves efficiently accelerate relativistic particles via 
resonant interaction. The condition for resonance between 
a wave of frequency u> and wavenumber projected along the 
magnetic field fcii , and a particle of type a with energy E a 
and projected velocity v\\ = vfi is (Melrose 1968; Eilek 1979): 



7 



- fciiwii = 



(15) 



where, in the quasi parallel case (k± « m a Q. a / p) , v — —1 
and v = 1 for electrons and protons respectively. 

The dispersion relation for Alfven waves in an isotropic 
plasma with both thermal and relativistic particles can be 
written as (Barnes & Scargle, 1973): 



l + ff72 + iic- 2 (l + jji) 



where 



m = 



^p,rcl 

nth 
and 



m e / n P;re i \ 



"2 = 



nth nip 



(16) 



(17) 



(18) 



Since the number density (and possibly the energy den- 
sity as well) of the thermal component in the ICM is con- 
siderably larger than the corresponding non-thermal com- 
ponent, one can show that r)i « 1 and r/2 « 1 so that the 
dispersion relation [Eq. (16)] becomes ui ~ |fe|||vA- Combin- 
ing the dispersion relation of the waves with the resonant 
condition, Eq. (15), one can derive the resonant wavenum- 
ber, k res , for a given momentum (p — mvj) and pitch angle 
cosine (u) of the particles: 

i, i fini 1 ,„„. 

k res ~ | = — -, (19) 



where the upper and lower signs refer to protons and elec- 
trons respectively. In an isotropic distribution of waves and 
particles, the particle diffusion coefficient in momentum 
space is given by (Eilek & Henriksen, 1984): 



D PP (p,t) 



2-K 2 e 2 v\ 
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W k (t) 



^ /t)A Amy 
V c pk J 



dk, (20) 



where the minimum wavenumber (maximum scale length) 
of the waves interacting with particles is given by : 



P 



1 

1± EA 



(21) 



is given by the largest wavenumber of the Alfven 
waves, limited by the fact that the frequency of the waves 
cannot exceed the proton cyclotron frequency, namely u) < 

fip. It follows that k m ax ~ Slp/VA Or k max ~ Qp/VM, V M 

being the magnetosonic velocity (here we assume k max ~ 

flp/VAl)- 

In the simple case of a power law spectrum of the MHD 
waves, W(k) oc k~ w (for k a < k < fc max ), Eq. (20) would 
give: 
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SB 2 = 8tt 



J W k dk, 
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(23) 



(24) 



w + 2 

The first order expansion (for va/c « 1) of Eq. 22 is 
the diffusion coefficient generally used in most recent the- 
oretical papers on electron acceleration in galaxy clusters 
(Ohno et al. 2002; Fujita et al. 2003). 



3.0.3 Electrons 

From Eq. (19), one can see that the momentum of the elec- 
trons which can resonate with waves with a given wavenum- 
ber k depends on the pitch angle cosine u. This resonant 
momentum can be written as 
Q e m e 1 



P : 



/' - 



(25) 



The minimum momentum of the electrons for which 
resonance with waves of a given wavenumber k can occur is: 



~k 



Q e + VAk 



which, in the relativistic limit becomes: 
rel _ n e m e 1 

Pmin i v A ' 



k l - 



(26) 



(27) 



Since the wavenumber of Alfven waves in a plasma is 
limited by ui < £2 P , from Eq. (26), one has that the minimum 
momentum of the electrons which can resonate with Alfven 
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va (m p 

Pmin =Pth — + 1 

Vth \ m e 



(28) 



which, in general, gives p min » p th , p th = m e v th being 
the momentum of the thermal electrons. It follows the well 
known result that thermal electrons cannot resonate with 
Alfven waves (Hamilton & Petrosian, 1992 and references 
therein). 

This important limitation of Alfven waves as particle 
accelerators forces us to consider the situation in which a 
relic population of relativistic electrons exists in the ICM 
(Sect. 2.2), and no electron acceleration from the thermal 
background is taken into account. 

In the case of an isotropic, homogeneus phase-space 
density electron plasma, the evolution of the electron spec- 
trum can be described in the context of the Fokker-Planck 
equation (Tsytovich 1966; Borovsky & Eilek 1986): 



df(p,t) = 
dt p 2 dp 



A> P P 2 ^ + f P 2 f(p) + 

Op dt rad 

dp 



+^.pf(p)\+ql(p,t), (29) 

where f(p, t) is the phase space density of electrons, D pp 
is the diffusion coefficient due to the interaction with the 
waves [Eq. (20)], dp/dt i and dp/dt rad give the ionization and 
radiative losses [Eq. (1-2)], and q£(p, t) is an isotropic phase- 
space electron source term. For simplicity we can introduce 
the two functions N(p,t) and Q e (p,t), related to / and q+ 
through the following relations: 

N(p,t)=4i:p 2 f(p,t), (30) 
and 

Q e (p,t) = 4irp 2 {qt(p,t)}. (31) 

The diffusion equation Eq. (29) in momentum space can 
therefore be transformed into an equation that describes the 
evolution of the electron number density: 



dN(p,t) _ d_ 
dt dp 



N(p,t) 



dp dp 

dt rad dt i 



d_ 

dp 



A: 



dN(p,t) 
dp 



+ Qe( P ,t). (32) 



3.0.4 Protons 

From Eq. (19) one has that the momentum of the protons 
which may resonate with waves with a given wavenumber k 
is: 



P - 



(33) 



The minimum momentum of the protons which may 
resonate with waves having wavenumber k is therefore: 



Pmin — (^p VAkJ . 

In the relativistic limit this reduces to: 



1 + 



(34) 



(35) 



From Eq. (34), one has that the minimum momentum of 
the protons which can resonate with Alfven waves in units 
of the momentum of the thermal protons is : 



va /fi p 

Pmin =Pth 

Vth V to 



(36) 



Since to < Q p , this basically means that thermal protons 
can efficiently resonate with Alfven waves (Hamilton & Pet- 
rosian, 1992). 

As in the case of relativistic electrons, the evolution of 
the proton spectrum is obtained from the equation: 



dN(p,t) _ d_ 
dp 



dt 



2 -D. 



d_ 

dp 



dN(p,t) 
dp 



+ Q P (p,t), 



(37) 



where dp/dt { is given by Eq. (5) and the coefficient D pp is 
provided in Eq. (20). 

3.1 From fluid turbulence to Alfven waves 

3.1.1 Injection 

We assume that fluid turbulence is present in the cluster 
volume with a power spectrum 



W t {xt) = W t °Xf 



(38) 



in the range x™ ln < Xf < a;f ax , where xf 1111 is the wavenum- 
ber corresponding to the maximum scale of injection of the 
turbulence and the maximum wavenumber is that at which 
the effect of fluid viscosity starts to be important and it is of 
the order of x?** ~ xf n (TZ)- 3/4 (e.g., Landau & Lifshitz, 
1959), 1Z being the Reynolds' number. 

For Kolmogorov turbulence we have m = 5/3, while 
for Kraichnan turbulence (Kraichnan 1965) one has m = 
3/2 (see Sect. 3.1.2 for a discussion on the Kolmogorov and 
Kraichnan phenomenology) . 

Here we investigate the connection between the fluid 
turbulence that we start with and the MHD waves that 
we use as particle accelerator. Fluid turbulence can radiate 
MHD modes (Kato 1968) via the Lighthill process. A fluid 
eddy may be thought of as radiating MHD waves in the 
mode j at a wavenumber k — (vf(x)/vj)x{, where Vj is the 
velocity of the j-mode wave. The MHD modes are expected 
to be driven only for 1 > it, it being the wavenumber at 
which the transition from large-scale ordered turbulence to 
small-scale disordered turbulence occurs. Following previous 
works in the literature (Eilek & Henriksen 1984; Fujita et 
al., 2003), we adopt the Taylor wavenumber as an estimate 
of this transition scale, namely: 



h = — 



2tt 

X T 



< ^ > / < b±r ) 



1/2 



dxi 
l (15/TZ) 1/2 , 



(39) 



where the Reynolds number is given by 1Z = l„»f/i*, and vk 
is the kinetic viscosity. The fraction of the fluid turbulence 
radiated is small for all but the larger eddies, near the Taylor 
scale, and thus the Lighthill radiation can be expected to not 
disrupt the fluid spectrum. More specifically, the energy rate 
radiated via the Lighthill mechanism into waves of mode 
j and wavenumber k is given by (e.g., Eilek & Henriksen 
1984): 



( k\ ~ y i 



(40) 
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where 



Vi = (2nj + 3) 



m — 1 



2n,(l - m) + 6 - 4m 



3 — m 



3 



(pu?Ji) 



(41) 



(42) 



and 



ft ~ pVf is the energy density of the fluid turbulence, 



R 



x Wf(x ) 
XTWf(xr) 



(43) 



Here, rij — for Alfven and slow magnetosonic (MS) 
waves, and rij = 1,2 for fast MS waves. In the following 
we concentrate on Alfven waves, in which case yj — 3/2 
(jjj — 1) for a Kolmogorov (Kraichnan) spectrum of the fluid 
turbulence. A more general treatment, including the effect 
of fast magnetosonic waves, will be given in a forthcoming 
paper. 

3.1.2 Basic equations and time evolution 

In our calculations we assume for simplicity that Alfven 
waves propagate isotropically in the cluster volume and we 
assume k ~ |fc|||. The spectrum of Alfven waves driven by 
the fluid turbulence evolves as a result of wave-wave and 
wave-particle coupling. In particular, the wave-particle in- 
volves the thermal and relativistic particles, in the way ex- 
plained in the previous Section. The combination of these 
processes produces a modified, time-dependent spectrum of 
Alfven waves, Ww(t), which can be calculated by solving the 
continuity equation (e.g., Eilek 1979): 



dW k (t) 
dt 



-&(D-^)-±*nm + w. («) 



The first term on the right hand describes the wave-wave 
interaction, with diffusion coefficient D kk = k 2 /r s . t s is the 
spectral energy transfer time, that for a given wavelength 
is given by r s ~ t% l /t3 (Zhou & Matthaeus 1990), where 
tatl = X/Sv is the non-linear eddy-turnover time (8v is the 
rms velocity fluctuation at A) and T3 is the time over which 
this fluctuation interacts with other fluctuations of similar 
size. 

In the framework of the Kolmogorov phenomenology, 
the Alfven crossing time ta = X/va largely exceeds tnl, 
therefore fluctuations of comparable size interact in one 
turnover time, namely T3 ~ tnl- In the Kraichnan phe- 
nomenology, ta << tnl, therefore convection limits the 
duration of an interaction and T3 ~ ta- Since the velocity 
fluctuation, 8v, is related to the rms wave field, SB, through 
the relation Sv 2 /v A = 8B 2 /B 2 , the diffusion coefficient is 
given by (Miller & Roberts 1995): 



£>kk 



t.7/2 / w k (*) \ 1/2 

,4 ( W k (t) \ 



(Kolmogorov) 
(Kraichnan) 



(45) 



in the Kolmogorov and Kraichnan phenomenology, respec- 
tively. 

The second term in Eq. (44) describes the damping with 
the relativistic and thermal particles in the ICM. In the 
case of nearly parallel wave propagation (k± << mfl/p, 
k ~ |fe|||) and isotropic distribution of particles of type a, 



the cyclotron damping rate for Alfven waves is as given by 
Melrose (1968): 



r£(t) 



47v 3 e 2 v 2 A 
kc 2 



„2„2 2 rPn 
7r e « A ' 

kc 2 



(1 - 



p 2 (l 

,N a (p,t) 



2 s df a (p, t) 

Ma) Tj~ dp - 

dN a (p,t) 



dp 



where, for relativistic particles, one has 



/'a 



VA 



pk 



while for sub-relativistic particles: 
th «Am . ric«m Q 



I'a 



V 



pk 



(46) 



(47) 



(48) 



Here the upper and lower signs are for negative and positive 
charged particles respectively. Lacombe (1977) showed that 
the damping rate for isotropic Alfven waves are well within 
a factor of ~ 3 of that calculated for nearly parallel wave 
propagation, therefore we are justified to use Eq. (46) in our 
calculations. 

The third term in Eq. (44) describes the continuous 
injection of Alfven waves as radiated by the fluid turbulence 
through the Lighthill mechanism. From Eq. (42) with rij = 
one has: 



h 



3- 2m 



3 — m 



3 



(jLy 



x k 



(49) 



4 GENERAL RESULTS ON THE EVOLUTION 
OF PARTICLES AND WAVES 

The interaction between waves and particles is investigated 
by solving the set of coupled differential equations given 
in Eqs. (32), (37), and (44). More specifically, the time- 
evolution of the wave spectrum depends on the damping pro- 
cesses which are affected by the spectra and energy content 
of the electron and proton components in the ICM. In turn, 
the wave-particle interactions modify the spectra of elec- 
trons and protons while energizing these particles. The var- 
ious components are clearly strongly coupled to each other 
and they cannot be treated separately as has been done in 
previous work on the subject. 

In the following we discuss separately in some detail the 
processes of turbulent cascading and damping. 



4.1 Turbulent cascade 

The equation that describes the turbulent cascading without 
accounting for damping processes and wave injection is the 
following: 



dw k (t) _ d t dw k (t) \ 
dt dk\ kk dk ) 



dw k (ty 



(50) 



The cascade timescale at a given wavelength is r kk 
fc 2 /D kk , and using Eq. (45): 
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Log[l(pc)] Log[l(pc)] 



Figure 5. Left Panel: Time evolution of the spectrum of Alfvcn waves injected in a single burst at a given scale. The spectra are 
plotted for 10 14 (solid line), 5 X 10 15 (dotted line), 3 X 10 16 (dashed line), and 3.2 X 10 16 s (dot-dashed line) after the injection event. 
In the calculations, a Kolmogorov diffusion coefficient is adopted. The temperature of the gas, the magnetic field and the gas density 
are T = 10 8 K, B = l^tG, and nth = 10~ 4 cm~ 3 respectively. Right Panel: Time evolution of the cascade— time scale estimated for the 
same parameters and with the corresponding line-styles as in the Left Panel. 



,,, 2 X 10 V fhoo 
Tkk(t) ^ 5 I : 

B^g V kpc 



)( 



nth 

10- 3 



' v^(^)" 1 
(Kolmogorov) 

(Kraichnan) 



4.2 Damping processes 

The first obvious damping process is provided by the inter- 
(51) action of the Alfven waves with the protons in the thermal 
plasma. Assuming that thermal particles have a Maxwellian 
energy distribution with temperature T and energy density 
£th, the damping rate is obtained from Eq. (46): 



where we define 5B>k ~ -\/8irkWk- It is worth noticing that 
the cascade timescale in the Kolmogorov regime does not 
depend on the value of the magnetic field strength. We also 
notice that in both Kolmogorov and Kraichnan regimes, the 
cascade timescale depends on the scale of the waves and 
on the density of the thermal plasma. In particular, Tkk is 
smaller in low density regions. 

In Fig. 5a we plot the evolution of the spectrum of a 
population of waves injected at a given scale as obtained 
from Eq. (50) in the Kolmogorov phenomenology (see cap- 
tion for details); the broadening of the wave distribution at 
scales larger than the injection scale clearly shows the effect 
of stochastic wave-wave diffusion. In qualitative agreement 
with Eq. (51), it is clear from Fig. 5a that the developing 
rate of wave-wave cascade from larger to smaller scales in- 
creases with decreasing scale. In Fig. 5b we show the cascade 
time-scale [from Eq. (51)] calculated for the spectra (and 
corresponding times) used in Fig. 5a. As a general remark, 
we find that for typical conditions of the ICM, the wave- 
wave time scale below 1 pc, namely on the scale relevant 
for wave-particle interaction, is considerably shorter than 
10 7 yr. 



r p+ ( fe ) = 

p th 

,2. .2 



lQ-K 3/2 e 2 vlm v £, 



3c 2 p? h fc 



1 + 



kVAJ 



m 



pVA ( 



k 2 v\ 



3c>? h fc 



exp 



exp 



V Pth / 

/ _ / Pmin s2\ 

V Pth / 

V Pth / 



where 



Pth = (2m p k B T) 



1/2 



(52) 



(53) 



and Pmin is given by Eqs.(34). 

It is well known that since the resonant condition [Eq. 
(15)] selects the interaction between particles of momentum 
p with waves with wavenumber k oc p~ , most of the en- 
ergy of the waves is dissipated in the acceleration of rela- 
tivistic particles (e.g., Eilek, 1979), and more specifically of 
relativistic protons. Although the exact damping rate with 
cosmic ray protons used in our calculations is obtained by 
numerically combining Eq. (46) and the solution of Eq. (37) , 
it is also possible to write an analytical expression for the 
damping rate for the asymptotic solutions obtained for the 
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proton spectrum [Eqs. (12-13)]. For supra-thermal protons, 
for /3 < Xm, the damping rate is given by: 



r P+ (fc) 



27r 2 e 2 m p K p v\ 



nthC 



C k s{s- 1) 



/n th (l- S )TCN| 
V x^mpC J J 



(54) 



while, for /3 > .X'n 



r P+ (fc) 



n t h (m p c 2 ) 2 C fc 
x 



r- 



2 (\l v m v 
(s-l)(s-3)V 



n p m 2 p v A r 1 



-(3-1)731 



s-3 



^-(.-^/s-ii 



3 dp}. (55) 



Here 



C = 1 + 



3CTn th 
p 3 /(m p c) 2 



(56) 



Finally, in the ultra-relativistic limit we distinguish the two 
regimes pi ow > p min (fc) and pi ow < p min (fc), where pi ow is the 
low energy cutoff in the proton distribution and p m \ n (k) is 
given by Eq. (35). 

For pi ow > Pmin(fe) the damping rate is given by : 



r P+ (fc) 



iv 2 e 2 v 2 A £ p (s - 2) 1 



k 



s + 2 



+ 



2(s + 2) v A &r,m 



s 

p i up 



S + 1 C Pl ow fc 

while, for pi ow < p m in(fe), one has 



(57) 



r P+ (fe) = 



27T 2 vi 



/ Plow 

\ /Tip 



Op 

IP 



(s + 2)(s-2) 



s + 1 



A;' 5 



(58) 



The damping rate due to relativistic electrons is ob- 
tained combining Eq. (46) and the solution of Eq. (32). A 
first comparison between the strengths of the wave-damping 
rates due to relativistic electrons and protons can be ob- 
tained assuming a simple power law energy distribution of 
the relativistic particles. For k « k m ax, the resonance oc- 
curs for pi ow < Pmin(fc) (i.e., Eq. (58) for protons and elec- 
trons) for both species and the ratio between the proton and 
electron damping rate is given by : 



F p + 



-(f) 

^low 



Plow 
Plow 



(59) 



where the signs + and — refer to the case of protons and 
electrons respectively. From Eq. (59) it is immediately clear 
that (for a typical p+ w /p{ ov , > 100, and for A^ 1 ~ A e rcl ) 
the damping rate on protons largely dominates that on elec- 
trons. 



1 2 



$0 




— 3 — 2 
Log(k/k max ) 



Figure 6. Comparison between damping and cascade time- 
scales. The plot shows the time scale for damping on the thermal 
gas (dashed line), on relativistic electrons (dotted line), and on 
relativistic protons (solid line). The cascade time scale is plot- 
ted as a solid straight line. Calculations are carried out assum- 
ing a Kolmogorov diffusion coefficient, and adopting T = 10 8 K, 
nth = 10~ 3 cm~ 3 , S e = 10~ 3 X £ th , £ p = 10~ 2 X £ th , s = 2.2 and 
Zi = 1.0, d(5B) 2 /dt = 3.3 x 10- 15 (u.G) 2 /s and I k oc fc- 3 / 2 . 



4.3 Damping versus Cascading 



The global damping time can be written as 

3 

".i 



- (£«)"'■ 



(60) 



For typical conditions in the ICM, the damping time on the 
thermal proton gas is < 10 5 sec (but the process is efficient 
only for fc/fc max > 0.1). The damping time on the relativistic 
component (especially protons) is usually > 10 8 sec. 

The time scale for the development of the wave-wave 
cascade depends on the wave- wave diffusion coefficient, Dkk, 
and thus on the energy density of the waves (Sect. 4.1). 
Given a spectrum of injection of waves per unit time, Ik, 
one simple possibility to estimate the cascade time scale 
and thus to compare it with the time scale of the damping 
processes is to use the spectum of the waves under stationary 
conditions and without damping processes, namely 



(ir-jf J ' (Kolmogorov) 



1/2 



(61) 



£r-£r) , (Kraichnan) 



The wave-wave time scale is therefore given by 

1/3 



k Z 1 

D kk k 



' (€) /(« /3 ), (Kolmogorov) 



. (€) 1/2 /(-Tl 1 k /2 ), (Kraichnan) 



(62) 
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A comparison between the time scales of the damping 
processes and of the wave- wave cascade is given in Fig. 6 for 
typical values of the parameters (assuming a Kolmogorov 
phenomenology, see caption). Fig. 6 shows that the time 
scale due to the damping with the thermal pool is consider- 
ably shorter than the cascade time scale for k/k max » 0.1 
so that a break or a cutoff in the spectrum of the waves is 
expected at large wavenumbers. However, the most impor- 
tant result illustrated in Fig. 6 is that, if a relatively large 
number of relativistic protons is present in the ICM, the re- 
sulting damping time scale can become comparable with or 
shorter than the wave-wave cascade time scale. This means 
that, at the corresponding wavenumbers, the spectrum of 
the waves is modified by the effect of the dampings and 
therefore that a power law approximation for the spectrum 
of the MHD waves cannot be achieved. We also note that 
the effect of the damping due to relativistic protons is par- 
ticularly evident at those wavenumbers that can exhibit a 
resonance with the bulk of the relativistic electrons in the 
ICM (those with 7 ~ 200 - 1000) and thus that this effect 
may have important consequences for the acceleration of the 
relativistic electrons. 

4-3.1 Dependence on the spectrum and energetics of 
protons 

As already mentioned in the previous Section, the damping 
of the Alfven waves on the relativistic protons modifies the 
spectrum of the waves and therefore indirectly affects the 
acceleration of electrons. We discuss here in some detail the 
dependence of this damping process upon the spectrum and 
energetics of the proton component in the ICM. 

The damping of the waves at a given wavenumber ba- 
sically depends on the number of protons with momentum 
that can resonate with such waves. At fixed number of rel- 
ativistic protons with supposedly a power law spectrum 
N(p) oc p~ s , the damping rate at wavenumbers correspond- 
ing to p » piow (piow being the minimum momentum in 
the proton spectrum) decreases with increasing s. 

The efficiency of the damping in galaxy clusters is fur- 
ther reduced in the case of steep spectra since, in this case, 
Coulomb losses affect the bulk of protons. 

The computed damping time scales, obtained using Eq. 
(60), are plotted in Fig. 7 for different proton spectra con- 
taining the same number of injected protons (see caption): 
only for large wavenumbers, where the condition p » pi ow 
is not applicable and the spectrum of the resonating pro- 
tons is not a power law, the damping rate in the case of 
steep spectra is comparable to that obtained for flat proton 
spectra. 

4-3.2 Dependence on the physical conditions in the ICM 

The development of the turbulent cascade and the damp- 
ing of the waves are affected mainly by the density of the 
thermal plasma n th , and by the magnetic field in the ICM. 

i) From Eq. (46) and Eqs. (54-58), at the zeroth order, 
one has that the damping time scale due to protons, Eq. 
(60), scales as Td oc £~ nth, while from Eq. (62) one has 
that the time scale for wave- wave cascade is r kk oc n\(^ . 
Both these time scales increase with increasing density of 




1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 

-5 -4 -3 -2 -1 
Log(k/k max ) 

Figure 7. Time scale for damping of Alfven waves on relativistic 
protons. From bottom to top the curves show this time scale for 
the following choices of the slope of the proton injection spectrum 
and energy content (e.g., Fig. 4) with respect to the thermal bud- 
get: s = 2.2 and £ p = 0.1 X £ th , s = 2.6 and £ p = 0.016 X £ th , 
s = 3.2 and £ p = 0.002 X £ th , s = 3.8 and £ p = 0.0005 X £ th . For 
comparison the cascade time-scale is plotted as a solid straight 
line. 



the ICM. However, the damping time due to protons is more 
sensitive to the gas density, in a way that, for a given £ p , the 
damping processes become less important for the shape of 
the spectrum of the waves as n t h increases: the comparison 
between the two time scales for two values of n t h is shown 
in Fig.8. 

If one assumes now that £ p oc n t h, the situation may 
change: in this case the damping time due to protons is 
constant and the importance of the damping processes is 
minimized in a low density ICM. 

ii) The effect of the magnetic field strength on the ef- 
ficiency of the damping processes and on the wave-wave 
cascade is more complex to evaluate, also due to the fact 
that the maximum wavenumber of the MHD waves de- 
pends on the magnetic field strength (in this paper it is 
kmax oc B). Assuming a time-independent injection power 
of MHD waves, J Ikdk, from Eq. (62) one has that the time 
scale for the wave-wave cascade is Tkk oc B -1 ' 3 . Such a sim- 
ple dependence is not obtained for the damping processes 
since from Eq. (46) and Eqs. (54- 58) one has that for large 
values of the wavenumber k, the damping time scale de- 
creases with increasing B, while for smaller values of k, 
it is Td oc B. This is illustrated in Fig. 9: the efficiency of 
the damping on ultra-relativistic protons (at small k) with 
respect to the turbulent cascade decreases with increasing 
magnetic field strength. On the other hand, the opposite 
trend is seen in the case of mildly and trans-relativistic 
protons (i.e., k/k max < 10~ 3 ). In general, the accelera- 
tion of relic (7 < 1000) relativistic electrons is powered by 
the resonance with waves with relatively large values of the 
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Figure 8. Time scale for damping of Alfven waves on protons as 
compared with the cascade time scale (straight lines) for nth. = 
10~ 4 cm -3 (solid lines) and 10 _3 cm -3 (dashed lines) and for a 
constant energy density in the form of cosmic ray protons. The 
injected spectrum of cosmic ray protons is taken as a power law 
with slope s = 3.0. The injection starts at z, = 1.0. The energy 
density is £ p = 0.2 X £ t h and = 0.02 X £ t h in the low and high 
density case, respectively. All the other parameters are fixed as 
in Fig. 6. 

wavenumber (e.g., k/k max > 1CP 3 ). Thus, given an injection 
rate of Alfven waves, the results shown in Fig. 9 indicate that 
the waves necessary to accelerate electrons with 7 < 10 J get 
appreciably damped by relativistic protons when the mag- 
netic field strength in the ICM is increased. The opposite 
trend is seen for the waves which resonate with very high 
energy electrons (e.g., 7 > 10 5 ). On the other hand, how- 
ever, it should be stressed that, in general, the efficiency 
of the wave-particle resonance increases with increasing B 
[Eq. (20, 22)], thus the conclusions given in this paragraph 
do not automatically imply that the electron acceleration is 
more efficient in regions with low magnetic field. 



5 QUASI STATIONARY SOLUTIONS 

For a given model for the injection of the waves in the ICM, 
the spectra of electrons, protons and waves can be calculated 
using Eqs. (32) with Q e — 0, (37) and (44). A few comments 
are in order: 

i) In this paper we confine our attention to the case of 
Alfven waves, that resonate predominantly with relativistic 
particles. Moreover, we do not consider situations in which 
the amount of energy injected in the form of turbulence 
becomes comparable with the thermal energy of the ICM. 
As a consequence, we can safely assume that the thermal 
distributions of electrons and protons in the ICM are not 
appreciably affected by the acceleration processes discussed 
above. 




B = 0.1mG 
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Figure 9. Time scale for damping of Alfven waves on relativistic 
protons compared with the cascade time scale for different values 
of the magnetic field strength: B = 0.1 (solid lines) and 2.0 fiG 
(dashed lines). It is assumed that n t ^ = 10~ 3 cm~ 3 , s = 3.0, 
£ p = 0.02 X Sth- Other physical parameters being are in Fig. 8. 
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Figure 10. Temporal evolution of the time scale of damping of 
Alfven waves on relativistic protons compared with the cascade 
time-scale and with the time scale for damping on the thermal 
protons (thick curves). The curves refer to for 0, 2 X 10 15 , 5 X 10 15 , 
1 x 10 6 sec, and 2 X 10 16 after the be ginning of the acceleration 
phase (from the top to the bottom (for fc/fc max < 0.1)). In the 
calculations we assume : d(SB) 2 /dt = 3.3 X 10" 15 (^G) 2 /s, T = 
10 8 K, n th = 10" 3 cm- 3 , B = 0.5/iG, £ e = 0.001 x £ th , £ p = 
0.0025 X £ t h, s = 3.0, pi n j > 0.1m p c, and Zj = 1.0. 



© 0000 RAS, MNRAS 000, 000-000 



Alfvenic reacceleration of relativistic particles in galaxy clusters: MHD waves, leptons and hadrons 13 




Figure 12. Acceleration time— scale (thin lines) and time— scale for energy losses (thick lines) for relativistic electrons as a function of 
the Lorentz factor. From bottom to top, the acceleration time— scales are calculated at 2 X 10 15 , 5 X 10 15 , 10 16 , and 2 X 10 16 sec after 
the beginning of the acceleration. The following values of the parameters are adopted: d(8B) 2 /dt = 3.3 X 10~ 15 (/iG) 2 /s, T = 10 8 K, 
n th = 10~ 3 cm- 3 , B = 0.5/jG, and £ e = 0.001 X £ th . Left Panel: £ p = 0.2 x£ th , s = 2.0 and z< = 1.0; Central Panel: £ p = 0.025 X £ th , 
s = 3.0 and z t = 1.0; Right Panel: S v = 0.002 X £ th , s = 4.0 and * 4 = 1.0; 




Figure 13. Acceleration time-scales (thin lines) and losses time-scales (thick lines) for relativistic electrons as a function of the Lorentz 
factor for 2 X 10 14 sec (Left Panel) and 10 16 sec (Right Panel) from the beginning of the acceleration. Calculations are shown for 
nth = 3 X 10 _3 cm — 3 (red lines) and n t ^ = 10 — 4 cm -3 (blue lines), and for B = 0.1/^G (solid lines) and B = 1.5/jG (dashed lines). In 
the calculations we assumed: d(SB) 2 /dt = 3.3 X 10- 15 (AtG) 2 /s, T = 10 8 K, £ e = 0.001 X £ th , £ v = 0.01 X £ th , and s = 2.2, not specified 
parameters being those in Fig. 12. 



ii) The spectra of electrons, protons and waves, as dis- 
cussed above, result from a coupling between all these com- 
ponents: the spectrum of the waves develops in time due 
to the turbulent cascade until damping becomes efficient 
and particle acceleration occurs. It is worth noticing that 
the time scales for the processes of damping and cascading 
are quite different from those related to particle losses and 
transport. While the wave spectrum develops over ~ 10 7 



sec, particle acceleration occurs on time scales of > 10 1 sec, 
and we are interested in following the particle evolution for 
a typical time of > 10 15 sec. 

The clear difference among these time scales suggests that 
we use a quasi stationary approach, in which it is assumed 
that at each time-step the spectrum of the waves approaches 
a stationary solution (obtained by solving Eq. (44) with 
dW/dt — 0) and that this solution changes with time due 
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Figure 11. Temporal evolution of the spectrum of the Alfven 
waves at times 2xl0 15 , 5 X 10 15 , 8 X 10 15 , 10 16 , and 1.5 X 10 16 sec 
after the beginning of the acceleration (from top to bottom). 
The calculations are carried out for a Kolmogorov spectrum 
of the fluid turbulence (i.e., yj = 3/2), a Kolmogorov diffu- 
sion coefficient d(SB) 2 /dt = 3.3 x 10- 15 (AtG) 2 /s, T = 10 8 K, 
n th = 10 _3 cm- 3 , B = 0.5^G, £ e = 0.001 xS th , £ v = 0.005 xS th , 
s = 3.2, Zj = 1.0 and pi„j > 0.1 m p c. The Taylor scale is at 

k ^ 10 fcmax- 
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Figure 14. Time evolution of the spectrum of cosmic ray protons 
as a function of p. From bottom to top, the curves are obtained 
at times 2 x 10 15 , 5 X 10 1B , 8 X 10 15 , 10 16 , and 1.5 X 10 16 sec 
after the beginning of acceleration. The initial proton spectrum 
is plotted as a dashed line. The values of the parameters are as 
in Fig. 11. 



to the evolution of the spectrum of the accelerated electrons 
and protons. In other words, at any time-step we solve a set 
of three differential equations: Eq. (44) (with dW/dt = 0), 
Eq. (37) and Eq. (32). Intermittent injection of turbulence 
in the ICM may occur on time-scales > 10 7 — 10 8 yrs which 
are much longer than the time-scales of damping and cas- 
cading, and thus the quasi-stationary approach discussed 
above remains applicable. 



5.1 The spectrum of Alfven Waves 

The shape of the spectrum of waves at any time is deter- 
mined by the damping of these waves, mainly on protons. 
The proton spectrum in turn changes because of accelera- 
tion, and backreacts upon the spectrum of waves: this im- 
plies that even for a time-independent rate of injection of 
waves, the strength of the damping rates and the spectrum 
of the MHD waves are expected to change with time. 

In Fig. 10 we plot an example of the evolution of the 
time scale for the damping due to relativistic protons: it is 
clear how the damping rate increases with time, as a conse- 
quence of the fact that most of the energy injected in MHD 
waves is channelled into relativistic protons. 

A relevant example of the time evolution of the spec- 
trum of waves is illustrated in Fig. 11: as expected, the 
energy associated with MHD waves which contribute to 
the acceleration of the bulk of the relativistic electrons 
(k/kmax ~ 10 _1 — 10~ 3 ) decreases with time. In addition to 
the general finding that the spectrum of the Alfven waves 
evolves with time, here we also point out that : 

a) the spectrum is not a simple power law: it has a different 
slope at different k and the curvature of the spectrum also 
changes with time; 

b) the spectrum has a low- A; cutoff due to the maximum 
injection scale, close to the Taylor scale; 

c) the spectrum has a high-fe cutoff generated by the 
damping with the thermal particles. 



5.2 Electron acceleration 

The process of electron acceleration via Alfven waves has 
been extensively investigated in the literature. Eilek & Hen- 
riksen (1984) showed that a self-similar solution for the spec- 
tra of electrons and waves can be found if these spectra are 
both required to be power laws in energy and wavenumber 
respectively. In this case, the slopes of the electron (5) and 
of the waves (ui) spectrum are related by 8 = 6 — w. 

As pointed out above, the assumption of power law 
spectra is usually not fulfilled. In fact, more often a sta- 
tionary solution is found in the form of a pile up spectrum 
N(p) ocp 2 cxp{-p/p c } (e.g., Borowsky & Eilek, 1986). 

In the general case considered here, the electron spec- 
trum may be even more complex due to the fact that the 
spectrum of the waves is not a power law and the whole evo- 
lution is time-dependent. Note that the initial stage of reac- 
celeration of relic relativistic electrons (i.e. 7 ~ 100 — 1000 
electrons) is mainly affected by the competition between 
Coulomb losses and acceleration due to the Alfven waves, 
while later stages, of further acceleration to the highest al- 
lowed energies, are limited by radiative losses. 
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In Fig. 12 we plot the time scale for electron acceler- 
ation, ~ p 2 /2D pp , and compare it with the time scale of 
the energy losses of electrons for different proton spectra in- 
jected in the ICM (see caption). For steep proton spectra, 
the electron acceleration is more efficient because less energy 
gets channelled into the proton component. In general, hard 
proton spectra make the acceleration of electrons to Lorentz 
factors 7 > 10 3 relatively difficult. 

If this is true in the initial stage of evolution of the 
system, it becomes increasingly less so at later times: after 
about ~ 0.5 — 0.7 Gyr, relativistic protons have accumu- 
lated enough of the waves energy that the damping of the 
waves becomes even more efficient and further acceleration 
of electrons is prevented. 

In Fig. 13 we plot the acceleration and cooling time 
scales of electrons as computed for different values of the 
density of the ICM and of the magnetic field strength. 

It is clear that at the beginning of the reacceleration 
phase (Fig. 13a) the electron acceleration is enhanced by in- 
creasing the magnetic field strength and by decreasing the 
density of the ICM. On the other hand, at later times in the 
reacceleration phase, the situation can be much more com- 
plicated. In particular, when the reacceleration efficiency is 
very high (for instance for large values of the magnetic field 
strength and for low values of rith), the energy stored in rel- 
ativistic protons after few reacceleration times may become 
sufficiently high to increase the damping rate of the Alfven 
waves and thus decrease the efficiency of electron accelera- 
tion (Fig. 13b): in this case, higher values of the magnetic 
field strength produce a lower acceleration efficiency com- 
pared with the case of low magnetic field. 

The continuous backreaction between waves and pro- 
tons creates a sort of wave-proton boiler that in a way is 
self-regulated. 

If the injection of fluid turbulence is intermittent on 
time scales of the order of the cooling time of electrons with 
Lorentz factors 7 ~ 10 3 — 10 4 , then the effect of the wave- 
proton boiler on the electron acceleration may be reduced. 
The reason for this is that for a given reacceleration rate, 
the accumulation of energy in the form of relativistic protons 
requires longer times and the electron acceleration remains 
efficient for ~ 1 Gyr. 

5.3 Proton acceleration 

Alfvenic acceleration of thermal and relativistic protons is 
extensively studied in the literature, in particular as applied 
to the case of solar flares (e.g., Miller, Guessoum, Ramaty 
1990; Miller & Roberts 1995). 

We consider here the extension of these calculations to 
the case of Alfvenic acceleration in clusters of galaxies. 

If the energy injected in Alfven waves is significantly 
larger than that stored by the relativistic protons at the 
beginning of the acceleration phase, then the spectrum of 
protons is expected to be considerably modified. Under these 
conditions, we illustrate, in Fig. 14, the evolution of the 
spectrum of the relativistic protons. The figure clearly shows 
that the spectrum flattens and develops a bump. 

The prominence of this bump increases with time as 
the energy absorbed by relativistic protons also increases. 
Moreover, the bump moves toward larger momenta of the 
particles during the acceleration time. 



The presence of a bump in the proton spectrum may 
be of some importance in the calculation of the spectrum of 
the high energy secondary electrons which are expected to 
be produced during hadronic collisions. We find that, under 
typical conditions in the ICM and assuming an energetics 
of the Alfven waves considerably larger than that of the 
initial proton population, relevant bumps are produced at 
energies in excess of 200 GeV which should be visible in the 
spectrum of secondary electrons at 7 > 10 5 . We also point 
out that proton acceleration is strongly reduced at some mo- 
mentum p m ax (Fig. 14) which corresponds to the momentum 
at which the resonance condition [Eq. (15)] is satisfied for 
wavenumber corresponding to the maximum injection scale 
of the Alfven waves [Eq. (39)]. A maximum energy of the 
injected secondary electrons is expected as well. A detailed 
discussion of the production of secondary electrons by reac- 
celerated protons and of their further acceleration will be 
presented in a forthcoming paper. 

A relevant point to make for the process of proton ac- 
celeration is to quantify the fraction of the energy injected 
in Alfven waves which is absorbed by the relativistic pro- 
tons and the time needed for relativistic protons to react 
to the injection of this energy. In Fig. 15, we compare the 
energy injected in waves with that converted to relativistic 
protons, for different values of the rates of energy injection 
in the form of Alfven waves. After some time, it can be seen 
that the energy of the protons asymptotically converges to 
the energy of the waves. It is also clear that there is a delay 
time between the injection of waves and the proton acceler- 
ation, as it is expected due to the finite acceleration time. 

5.4 The Wave-Proton Boiler 

One of the most important results of our investigation is the 
quantitative treatment of the backreaction of the accelerated 
protons on the waves and in turn on the electrons. Qualita- 
tively, given typical conditions in the ICM, we can identify 
three main temporal stages of the acceleration process: 

1) Cascading stage: 

For a non negligible rate of energy injection in the form 
of Alfven waves, the cascade time is shorter than the damp- 
ing time. This remains true up to some critical wavenum- 
ber, which depends on energetics and spectrum of protons, 
where damping starts to be relevant. If such a wavenumber 
is larger than about W~ 2 k max , then enough energy is left 
in the form of waves at the scales which may resonate with 
relic relativistic electrons. In this case electrons are effec- 
tively re-energized. 

2) Stage of proton backreaction: 

Once the Alfven waves start to accelerate electrons and 
protons to higher energies, the spectrum of protons and elec- 
trons becomes harder and the fraction of the energy stored 
in non-thermal particles starts to be large enough to make 
damping more severe. As a consequence, the rate of electron 
acceleration is reduced. 

3) End of acceleration: 

At the beginning of the acceleration phase, the bulk of 
protons is located at supra-thermal or trans-relativistic en- 
ergies. It takes a few 10 8 yrs, however, for these protons to be 
energized to higher energies, as illustrated in Fig. 16, where 
we plot the acceleration time scale of relativistic protons 
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Figure 15. Evolution of the energy density stored in the form 
of cosmic ray protons (p > 0.1 m p c) as a function of time from 
the beginning of the acceleration stage for three different injec- 
tion rates of Alfven waves: d(6B) 2 /dt = 3.3 X 10~ 15 (points-solid 
line), 1.9 X 10~ 15 (points-dotted line), and 8.2 X lCr 16 (^G) 2 /s 
(points-dashed line). For comparison the energy density injected 
in Alfven waves is aslo plotted. In the calculations wc assume 
the following values of the parameters: B = 0.5/^G, T = 10 8 K, 
nth = 10~ 3 cm- 3 , £ e = 0.001 x £ th , £ p = 0.01 x £ th , s = 2.2 and 
Zi = 1.0. 

(we chose a relatively steep injection spectrum). From Fig. 
16 we see that after about 0.5 — 0.7 Gyr the acceleration time 
scale has increased by about one order of magnitude. At this 
point the acceleration stage of protons and electrons can be 
considered as concluded, unless the injection of turbulence 
occurs intermittently (see Sect. 5. 2). 

After the end of the third stage, the electrons cool due to 
radiative and Coulomb losses, while the Alfven acceleration 
is only able to prevent the thermalization of these particles 
maintaining their Lorentz factor around 7 ~ 100 — 1000. 



6 NON-THERMAL EMISSION FROM 
GALAXY CLUSTERS 

6.1 Cluster mergers and turbulence 

Major mergers are among the most energetic events in the 
Universe. Cluster mergers involve a collision of at least two 
subclusters with relative velocity of ~ 1000 — 2000 km s _1 . 
During these events, a gravitational energy in excess of 10 63 
erg is released through the formation of shock waves. Numer- 
ical simulations show that cluster mergers can generate rela- 
tively strong turbulence in the ICM (Norman & Bryan 1999; 
Roettiger et al., 1999; Ricker & Sarazin 2001). The bulk of 
the turbulence is most likely injected on scales > 100 kpc 
during the motion of the subclusters. Afterwards this turbu- 
lence eventually cascades toward smaller scales. As discussed 
in Sect. 3.1, when the turbulent cascade reaches scales close 




Log(p) 

Figure 16. Evolution of the proton acceleration time-scale dur- 
ing the acceleration phase as a function of p. From bottom to top 
the curves refer to 2 X 10 15 , 5 X 10 15 , 10 16 , and 2 X 10 16 sec after 
the beginning of the acceleration stage. The following values of the 
parameters have been adopted: d(8B) 2 /dt = 3.3 X 10~ 15 (/iG) 2 /2, 
B = 0.5/iG, T = 10 8 K, nth = 10- 3 cm- 3 , £ e = 0.001 x £ th , 
£ p = 0.0025 x £ th and s = 3.0. 

to the Taylor scale, a fraction of the energy flux of the fluid 
turbulence can be transferred to MHD waves which in turn 
can accelerate fast particles. 

For simplicity, we assume here that the bulk of the fluid 
turbulence in a given point of the cluster volume is injected 
at the scale l a , for a time n, of the order of the time necessary 
for the subclump to cross the scale l a : 

n(Gyr)~0.3e(^)(^p, (63) 

where v c is the velocity of the subclump in the host cluster 
and ^ is a parameter of the order of a few. Within these 
assumptions, the injection rate of energy in the form of fluid 
turbulence is given by : 

2.3 x IP' 27 / na \ / T \ / l Q \ - 1 / v c \ £ t 

Ff 1 lio^XioO UooJ \W)s^' (64) 

where £ th is the local energy density of the ICM in the form 
of thermal gas and £ t is that in the form of turbulence. The 
bulk of the fluid turbulence at the scale l then cascades to- 
ward smaller scales producing a spectrum of the fluid turbu- 
lence that we write as Wf(x) oc x~ m (Sect. 3.1). Assuming 
a Kolmogorov phenomenology for the wave-wave diffusion 
in fc— space, the time scale for the cascade can be estimated 
from Eq. (62) with I Xo ~ x~ l Ff. 

r s (Gyr)~0.2( — j (— j (— j * (-^-J .(65) 

In our simple approach, this is the time delay between the 
merger event and the development of the turbulence at small 
scales (and thus the production of MHD waves) . In the case 
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of Alfven waves in the framework of a Kolmogorov phe- 
nomenology, the power injected in waves is [Eq. (49)]: 



Pa 



I 



I(k)dk : 



1.5 x 1Q- 2 V vt V D -x^ nth W la y 1 

ni6 k liooj MlO^J 1 300 J (66) 

All the quantities involved in the calculation of the power 
radiated in the form of Alfven waves can be relatively well 
modelled. The only parameter which is very difficult to esti- 
mate is the value of the Reynolds number, TZ = l Vf /vk, due 
to the large uncertainties in the value of the kinetic viscos- 
ity vk = u p X e ff/3 (u p is the thermal velocity and A e // the 
effective mean free path of protons). For transverse drift of 
protons in a magnetic field B, the mean free path is given by 
X e ff ~ Ag/A c (e.g., Spitzer 1962) with A 9 and A c being the 
proton gyroradius and the mean free path due to Coulomb 
collisions respectively. In this case the kinetic viscosity is 
given by (Fujita et al. 2003, and references therein): 



2 



(67) 



The resulting Reynolds number is: 

which is extremely large and, most likely, should be con- 
sidered as an upper limit, due to the assumption of trans- 
verse drift of protons on the magnetic field lines. In gen- 
eral, diffusion of the protons along the magnetic field lines 
can substantially increase the value of X e ff and thus re- 
duce TZ. The Reynolds number has been roughly estimated 
in a number of astrophysical situations, being ~ 10 7 in the 
solar wind (e.g., Grappin et al 1982), ~ 10 11 in the extra- 
galactic radio jets imaged with the VLA (e.g., Henriksen, 
Bridle, Chan 1982), ~ 10 14 in the solar corona (e.g., Ofman 
& Aschwanden, 2002) and ~ 10 for the hot phase of the 
local interstellar medium (e.g., Armstrong, Rickett, Cordes 
1981). Following this last estimate, in our modelling we take 
TZ ~ 10 16 but also stress that, due to the poor dependence 
of Pa on 1Z, an uncertainty in TZ by six orders of magnitude 
implies only one order of magnitude change in Pa- 

In our simple model for the merger, it is easy to show 
that the condition Pa « Ff is satisfied, therefore the ap- 
plication of the Lightill scheme in our calculations appears 
to be self-consistent. 

The power radiated in the form of Alfven waves strongly 
depends on the velocity of the eddies in the fluid turbulence. 
In addition, assuming that the energy of the fluid turbulence 
is a fraction of the local thermal energy (i.e., £ t ~ pv^, and 
v/=const) and that the largest scale of the spectrum of tur- 
bulence, l , does not depend on the position in the cluster 
volume, we notice that the injected power in Alfven waves 
increases with increasing the number density of the ICM and 
with decreasing the strength of the local magnetic field. For 
instance, assuming a typical scaling law for the magnetic 



the power injected in the form of Alfven waves is expected 
to be slightly larger in the central regions of the cluster. 

A second important quantity in our calculations is the 
largest scale of the spectrum of the Alfven waves, kr ~ 



XTVf(xT)/vA- It is convenient to parametrize this quantity 
in terms of the wavenumber corresponding to the minimum 
scale of the MHD waves, k max ~ 0. p /vm, namely: 



(69) 



fcr _ TZ 16 3 / n th \ 2 /300\ 2 / v f \ / T_\ ? 
fc max ~ 2.4- 10U 10-s J Wo / MG V400Al08j ' 

Given the resonance condition [Eqs.(27) and (35)], the 
presence of a maximum scale in the wave spectrum, It ~ 
2ir/kT, implies a limit for the energy of the particles that 
can be efficiently accelerated: 



^(TeV)<24x(^)(^ 



1/2 



oc TZ 



-1/3 



- 1 / 2 , D 2 -1 



(70) 



This limit provides a relatively good estimate in the case 
of relativistic protons, which are basically loss-free particles, 
while the maximum energy of the electrons is driven by the 
competition between acceleration and loss terms. 

6.2 Constraining the model parameters 

In this Section we derive some constraints on the physical 
conditions in the ICM in order to obtain the reacceleration 
efficiency necessary to allow the production of the observed 
non-thermal emission. 

From Eq. (66) it is clear that a key parameter in the 
calculation of the efficiency of particle acceleration in our 
approach is given by : 



V400/ UO 16 



\30oy ' 



(71) 



which is very sensitive to the velocity of the eddies in the 
fluid turbulence. A constraint on this velocity can be ob- 
tained by requiring that the energy injected in the form of 
fluid turbulence {Et ~ x W{x ) ~ m p nthVf) does not exceed 
the thermal energy. From Eqs. (63) and (64) we obtain: 



1/2 



(72) 



An additional constraint to the parameter ^ comes from 
requiring that the energy stored in the form of relativistic 
protons is small enough to allow for efficient electron accel- 
eration. 

As pointed out above, after an acceleration time of the 
order of 0.5-0.7 Gyr, relativistic protons get a large fraction 
of the energy previously in the form of Alfven waves. We 
limit ourselves to cases in which Pa ■ n < 0.1 £ t h- From 
Eqs. (63) and (66) one has: 



•*T«-(^r'(£)&r- 



(73) 



Finally, we are interested in the production of relatively 
long- living (e.g., for > 0.3 Gyr) non-thermal phenomena, in 
order to have a chance of observing them in some clusters. 
From Eq. (63) we can set a limit on the parameter £: 



field in the cluster B oc n^ 3 , it would be Pa oc n^ 6 , namely £ > 



,(— v 

1000 J 1.300/ 



(74) 



The maximum energy, jmax, of the accelerated electrons 
is obtained by balancing energy losses and energy gains. 
In Fig. 17a we plot ^) m ax as a function of for different 
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Figure 17. Maximum energy attainable for relativistic electrons during the acceleration (only systematic terms are considered here) 
plotted as a function of *Sf. Left Panel: 7 max at 3 X 10 15 (dashed line) and 10 16 sec (solid line) after the beginning of the acceleration 
stage. The values of the parameters are: B = 1/iG, n t ^ = 10 _3 cm -3 , £ e = 0.001 X £ t /j, £ p = 0.003 X £ t f l , s = 2.2, and Zj = 1.0. Right 
Panel: 7max at 10 16 sec from the beginning of the acceleration stage for £ p = 0.003 X Sth (solid line) and £ p = 0.03 X £th (dashed line). 
Other parameters are as in the Left Panel case. 



acceleration times, AT, for a given set of values of the pa- 
rameters denning the environment of cluster cores. In Fig. 
17b we plot the same quantity for different values of the 
energy density in the form of relativistic protons, £ p . The 
injection spectrum of protons is taken as a power law with 
slope 2.2. 

In order to obtain j m ax » 1000, needed to explain 
the synchrotron emission at GHz frequency as well as the 
IC hard X-ray photons, we are forced to require that * > 1 
and £ p < 3% Sth- 

Such a stringent limit is the consequence of the effec- 
tive damping of Alfven waves upon the relativistic proton 
component, which inhibits the acceleration of electrons. 

In the periphery of the cluster, where the magnetic field 
is expected to be lower, the conditions to obtain high energy 
electrons are less stringent. It remains true however that 
no more than a few percent of the thermal energy of the 
cluster can be in the form of relativistic protons if we want to 
interpret the observed non-thermal phenomena as the result 
of radiative processes of high energy electrons accelerated 
via Alfven waves. 

A steeper injection spectrum of protons, containing the 
same energy density, makes the constraints found above even 
more stringent; however this would imply a very large energy 
injected in relativistic protons in the ICM. 

On the other hand, at given proton number density, a 
steeper spectrum contains less relativistic particles, which 
would allow for more efficient electron acceleration. 

Clearly, a crucial parameter in the modeling of the non- 
thermal phenomena in galaxy clusters is the strength of the 
magnetic field in the ICM. There is still debate on whether 
this field is of several jj,G or rather fractions of /iG: on one 
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Figure 18. Synchrotron cut-off frequency as a function of the 
magnetic field strength in the case of n th = 10~ 3 cm -3 (empty 
symbols) and n t ^ = 2 X 10~ 4 cra~ 3 (filled symbols). The results 
are reported after 7 X 10 15 sec of reacceleration. In the calcu- 
lations the following values of the parameters have been used: 
d{5B) 2 /dt = 3.3 x 10- 15 (/iG) 2 /s, T = 10 S K, £ e = 0.001 x £ th , 
£ p = 0.01 X £ th , s = 2.2 and z t = 1.0. 
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hand, if the hard X-ray excess is interpreted as a result of IC 
emission of relativistic electrons, the required volume aver- 
aged field is < 0.2-0.4n:G (e.g. Fusco-Femiano et al. 2002). 
On the other hand, Faraday rotation measurements (RM) 
of the radiation coming from cluster radio sources seem to 
require a magnetic field strength of ~ 4 — 8/iG (e.g., Clarke 
et al. 2001). A number of possibilities to reconcile these two 
predictions have been proposed in the literature, based on 
either choosing more realistic spectra of the radiating par- 
ticles or introducing a spatial distribution of non-thermal 
particles and magnetic fields (Goldshmidt & Rephaeli 1993; 
Brunetti et al. 2001a; Petrosian 2001; Kuo et al. 2003). 

It is also worth recalling that Faraday RM do not pro- 
vide a direct measurement of the magnetic field, but rather 
an estimate of a quantity which is a (not necessarily triv- 
ial, e.g. Clarke 2002) convolution of the component of the 
magnetic field parallel to the line of sight, of the electron 
density and that accounts for the topology of the magnetic 
field: Newman et al. (2002) showed how the assumption of 
a single-scale magnetic field leads to an overestimate of the 
field strength extracted from RM. Other authors have fur- 
ther discussed the influence of the power spectrum of the 
magnetic field topology in the ICM on the RM (Ensslin & 
Vogt 2003; Vogt & Ensslin 2003; Govoni et al. 2003). 

In order to illustrate the effect of these uncertainties on 
the conclusions inferred from our calculations, we evaluated 
the so-called synchrotron cut-off frequency for a cluster with 
temperature T = 10 8 K and with a relic electron and proton 
energy densities chosen as £ e = 0.001 x £ t h, £ P = 0.01 x £ t h 
(protons have a spectrum with power index s = 2.2; see cap- 
tion of Fig. 18 for additional information). The calculation 
is carried out for a dense region (nth = 10 _3 cm~ 3 [empty 
symbols]) and for a low density region (nth = 2 x 10 _4 cm~ 3 
[filled symbols]). Given the shape of the spectrum of the ac- 
celerated electrons, a synchrotron cut-off at > 300 MHz is 
required to account for the synchrotron radiation observed 
in the form of radio halos. 

Our conclusion, based on Fig. 18, is that in high density 
regions (n th ~ 10~ 3 cm~ 3 ) Alfvenic reacceleration of relic 
electrons cannot be an efficient process for B 3> 4/nG and 
for B <C 0.5piG. These constraints become less stringent in 
the case of low density regions. 

6.3 A simplified models for Radio Halos and Hard 
X— ray emission 

The best evidence for the diffuse non-thermal activity in 
clusters of galaxies is provided by the extended synchrotron 
radio emission observed in about ~ 30 massive clusters of 
galaxies (e.g., Feretti, 2002). A recent additional evidence 
supporting the existence of relativistic electrons is given by 
the hard X-ray tails in excess to the thermal emission discov- 
ered by BeppoSAX and RXTE in the case of a few galaxy 
clusters (e.g., Fusco-Femiano et al. 2002). The possibility 
that hard X-ray tails are due to IC scattering of the CMB 
photons is intriguing as, in this case, radio and HXR radi- 
ations would be emitted by roughly the same electron pop- 
ulation and thus the combination of radio and HXR data 
would allow us to infer an estimate of the volume- averaged 
magnetic field strength and of the energy density of relativis- 
tic electrons. Additional pieces of evidence for non-thermal 
phenomena are the so called radio Relics and the EUV ex- 
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Figure 19. Temporal evolution of the accelerated electron spec- 
tra after 0, 10 14 , 5xl0 15 , 7xl0 15 , 10 16 , and 1.2 X 10 16 scc from 
the beginning of the acceleration stage. The following values of 
the parameters have been used: d(SB) 2 /dt = 1.6 x 10~ 15 (/^G) 2 /s, 
B = 0.5/iG, T = 10 8 K, n th = ltr 3 cm- 3 , £ e = 0.001 x £ th , 
£ p = 0.01 X £ th , s = 2.2 and z t = 1.0. 

cesses whose origin may however be not directly connected 
to that of radio halos and HXR. Therefore these phenom- 
ena will not be modelled in the present paper. For a recent 
review on these arguments the reader is referred to Ensslin 
(2002, and ref. therein) and Bowyer (2002, and ref. therein). 

In this section we apply the formalism described in pre- 
vious sections in order to show that for the conditions re- 
alized in the ICM, Alfvenic reacceleration of relic electrons 
may generate the observed radiation, provided the energy 
content in the form of relativistic protons is not too large. 
For simplicity we neglect here the production and reacceler- 
ation of secondary products of proton interactions. 
Our simple model for the ICM assumes a /3-model (Cavaliere 
& Fusco-Femiano, 1976) for the radial density profile of the 
thermal gas in the ICM, in the form 

n th (r)=n th (r = 0)(l + (-f) , (75) 

where r c is the core radius and we adopted (3 — 0.8. The 
magnetic field is assumed to scale with density according 
with flux conservation: 

Based on the constraints given in Sect. 6.2 we adopt B(r = 
0) ~ 0.5 - 4/xG. 

Finally, we assume that the ratio between the energy 
density of the relic relativistic particles (at the beginning of 
the acceleration phase) and that of the thermal plasma is 
constant with distance throughout the cluster volume: 

£p[e] = £thVp[e], (77) 

where r] p [ e ] is a free parameter (n < 0.1). 
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Figure 21. Temporal evolution of the synchrotron (Left Panel) and ICS (Right Panel) spectra from the cluster model given in 
Sect. 6.3 (integrating the emissivities up to 5 X r c ). Spectra are shown at 5 X 10 15 , 1.2 X 10 16 and 1.7 X 10 16 sec from the beginning of 
the acceleration (from bottom to top). In the calculations £ e = 5 X 10 — 5 £th nas been adopted. Radio data are taken from Thierbach et 
al. (2003), EUV from Bowyer et al. (1999; here reported as an upper limit, see text), XHR from Fusco-Femiano et al. (2004). 
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Figure 20. Spectra of electrons accelerated for 1.2 X 10 16 sec at 
different distances from the cluster center: r =0.3, 0.6, 0.9, 1.2, 
1.4, 2.1r c (from top of the diagram). The central values assumed 
in the calculations are : ra t ^(0) = 1.5 X 10~ 3 cm~ 3 , B(0) = 1.5^iG, 
d(SB(0)) 2 /dt = 2.2 X 10" 15 (mG) 2 /s with the scaling laws given 
in Sect. 6.3, r c = 400 kpc. 



For simplicity we also assume that the maximum in- 
jection scale of the turbulence, the Reynolds number and 
the velocity of the turbulent eddies are independent of the 
location within the cluster volume. 



Using the scaling relationship for the magnetic field, Eq. 
(76), and the expression for the injection power in the form 
of Alfven waves, Eq. (66), we obtain: 



5/6 



(78) 



The spectrum of the relativistic electrons in the core 
region is plotted in Fig. 19. We can see that the bulk of 
relativistic electrons, initially at 7 ~ 10 2 , can be energized 
up to 7 ~ 10 4 for a relatively long time. 

Eq. (78) indicates that, in our simple approach, the 
power injected in the form of Alfven waves decreases with 
increasing distance from the cluster center. Since radio ha- 
los have a considerable size, it is needed to check that our 
model provides enough energy in the outskirts of clusters. 

In Fig. 20 we plot the electron spectra at different dis- 
tances from the cluster center (see caption). What emerges 
from this figure is that at large distances the effect of ac- 
celeration is even stronger than in the central region and 
the electron spectra peak at slightly higher energies than in 
the core. This is due to the fact that in the outskirts the 
damping rate is reduced more than the rate of injection of 
turbulence. 

The synchrotron emissivity roughly scales as J syn oc 
N e j m ax(r) 2 B 2 oc n 7 / h 3 'j ma x(r) 2 . Such a relatively soft de- 
pendence of the emissivity on the radius allows for a rela- 
tively broad synchrotron brightness profile. 

More specifically, the synchrotron emissivity can be 
written as: 



Jsynij^i i) 



V3e 3 B 



0N(p,t)F{£), (79) 
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where the synchrotron Kernel is given by (e.g., Rybicki 
& Lightman, 1979): 



'(-) = - / 



(80) 



and v c — (3/47r)p 2 e_B sin 0/(mc)' i . 

The emissivity due to inverse Compton scattering off 
the CMB photons is given by (Blumenthal & Gould, 1970): 



~(yi,t) = (2nr m 2 c) hv\ I I dvdp 

J v J p 



N(p,t)p- 



Ap 2 v 



2p 



where T z = 2.73(1 + z) is the temperature of the CMB 
photons. 

The corresponding synchrotron and IC spectra are plot- 
ted in Fig. 21 at different times, for a central magnetic field 
B(r = 0) ~ lfiG. The spectra are compared with that ob- 
served for the radio halo in the Coma cluster : an initial 
energy density in the relic relativistic electrons of the order 
of 5 x 10~ 5 £ t h is required to account for the data. In Fig. 21 
we also report the luminosity of the EUV excess in the Coma 
cluster as an upper limit since it is commonly accepted that 
the origin of the EUV excess is not directly related to the 
same electron population responsible for the radio and pos- 
sibly for the HXR emission (Bowyer & Berghofer 1998; En- 
sslin et al. 1999; Atoyan & Volk 2000; Brunetti et al. 2001b; 
Tsay et al. 2002). Here we stress that Fig. 21 does not show 
the best fit to the data but just a comparison between data 
and time evolution of the emitted spectra resulting from the 
very simple scaling of the parameters described above. On 
the other hand, it should also be stressed that the time evo- 
lution of the synchrotron and IC spectra reported in Fig. 
21 is generated via the first fully self-consistent calculation 
of particle acceleration in galaxy clusters. Thus provided 
that the energy of relativistic protons in galaxy clusters is 
not larger than a few percent of the thermal energy, Fig. 
21 proves, the possibility to obtain the observed magnitude 
of the non-thermal emission in these objects via Alfvenic 
acceleration. 

In passing we note that the integrated radio spectrum 
at ~ GHz frequencies is steeper than the IC spectrum in the 
hard X-ray band independently of the time slice. This seems 
in nice agreement with some temptative evidences published 
in some recent literature (Fusco-Femiano et al., 2000; 2002). 



7 CONCLUSIONS 

The origin of the extended non-thermal emission in galaxy 
clusters is still a subject of active investigation. Despite the 
numerous models put forward to explain the observed fea- 
tures of this diffuse emission, only some of them appear to 
be able to describe the observations in detail. One of the 
approaches that has definitely produced the most promising 
results consists in the continuous reacceleration of relic rela- 
tivistic electrons leftover of the past activity occurred within 
the ICM. The reacceleration is likely to occur through reso- 
nant interaction of electrons and MHD waves. 



In this paper we presented a full account of the time- 
dependent injection of fluid turbulence, its cascade to 
smaller scales, the radiation of Alfven waves through the 
Lighthill mechanism, the resonant interaction of these waves 
with electrons and protons (namely their acceleration) and 
the backreaction of the accelerated particles on the waves. 

The solution of the coupled evolution equations for the 
electrons, protons and waves revealed several new interest- 
ing effects resulting from the interaction among all these 
components: 

i) Alfven waves are radiated by the fluid turbulence 
through the Lighthill mechanism. This allows us to establish 
a direct connection between the fluid turbulence likely to be 
excited during cluster mergers, and the MHD turbulence 
that may resonate with particles in the ICM. 

ii) Previous calculations looked for self-similar solutions 
for the spectra of electrons and MHD waves in the form 
of power laws. The solutions obtained in the present pa- 
per show that in Nature these self-similar solutions are not 
necessarily achieved. In general the system evolves toward 
complex spectra of electrons and MHD waves with a bump 
in the electron spectrum, that moves in time toward an in- 
creasingly large particle momentum. 

Hi) We limit ourselves here with the case of Alfven 
waves, which couple efficiently with both thermal and rel- 
ativistic protons, but only with relativistic electrons. The 
main damping of Alfven waves occurs on relativistic pro- 
tons, if there are enough of them. The spectrum of the waves 
is cutoff at small scales due to the damping of these waves. 
The damping moves energy from the waves to the particles, 
determining their acceleration/heating. 

iv) The importance of the presence of the relativistic 
protons for the acceleration of electrons is one of the most 
relevant new results of this work. A large fraction of the ther- 
mal energy in the form of relativistic protons enhances the 
damping rates of Alfven waves, suppressing the possibility 
of resonant interaction of these waves with electrons. Since 
electrons are the particles that radiate the most, a too large 
fraction of relativistic protons suppresses non-thermal phe- 
nomena directly related to electron reacceleration via Alfven 
resonance. It is worth reminding that in principle the sec- 
ondary electrons and positrons resulting from hadronic in- 
teractions may also be re-energized by waves. We neglect 
this effect here. Our results show that no more than a few 
percent of the thermal energy density of the cluster can be 
in the form of relativistic protons if we want to interpret 
the diffuse radio and hard X-ray emissions as the result of 
synchrotron and ICS radiation of relic electrons reacceler- 
ated through Alfven waves. This appears as a stringent con- 
straint on the combination of proton number and spectrum 
since the accumulation of large number of protons in the 
ICM is predicted by both analytical calculations (Berezin- 
sky, Blasi & Ptuskin 1997) and numerical simulations (Ryu 
et al. 2003). The amount of energy piled up in the form of 
protons in the ICM depends strongly upon the history of 
cosmic ray injection: the ICM is polluted with cosmic ray 
protons both due to single sources such as galaxies, radio 
galaxies and active galaxies, but also due to diffuse accel- 
eration events, such as mergers of clusters of galaxies and 
accretion of cosmological gas onto a gravitational potential 
well of a cluster which has already been formed. The spec- 
trum of the cosmic rays diffusively trapped in the ICM de- 
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pends upon their origin. It is expected that mergers may 
contribute a major fraction of these cosmic rays, but the 
spectrum has been calculated to be relatively steep because 
of the weakness of the merger related shocks. Accretion on 
the other hand, contributes spectra which are as flat as E~ 2 . 
Both the energy and spectrum of the cosmic rays stored in 
the ICM affect the temporal evolution of Alfven waves and 
their ability to reaccelerate relic electrons. If future observa- 
tions will unveil the presence of a population of relativistic 
protons in the ICM with > 5 — 10% of the thermal en- 
ergy, then Alfvenic reacceleration of relativistic electrons in 
galaxy clusters will be discarded as a possible explanation 
of non-thermal phenomena in the ICM. On the other hand, 
similar reacceleration phenomena can be driven by MHD 
waves other than Alfven waves (e.g., MS, LH), for which 
the energy transfer does not occur preferentially toward pro- 
tons. In this case, the bounds presented here on the allowed 
energy density in the form of relativistic protons in the ICM 
could be substantially relaxed. 

v) By assuming that protons have a relatively flat spec- 
trum (s ~ 2.2) and that they contain up to a few percent 
of the thermal energy, we used a simple but phenomeno- 
logically well motivated model for the density and magnetic 
field in a cluster of galaxies in order to calculate the expected 
non-thermal radio and hard X-ray activity of the cluster as a 
function of time. For the first time we performed a fully self 
consistent calculation and showed that the reacceleration 
of relic electrons through resonant interaction with Alfven 
waves can explain very well the observed phenomena, in- 
cluding the extended diffuse appearance of this emission. In 
passing, without considering the case of the HXR emission, 
we also showed that a magnetic field strength in the range 
0.5 — 4/iG in the cluster cores allow electron acceleration ef- 
ficient enough to produce GHz synchrotron emission. These 
conditions are less stringent in the outermost regions of the 
clusters. 

vi) In our calculations we adopted a constant injection 
rate of turbulence during particle acceleration. In this case, 
even assuming £ p < 0.1 x £ t h, we find that after a few 10 8 yrs 
of acceleration the backreaction of protons on MHD waves 
can suppress the acceleration of energetic electrons : this 
provides a limit on the duration of the non-thermal phenom- 
ena in galaxy clusters. It is possible to extend the duration of 
non-thermal activity assuming that injection of turbulence 
occurs in relatively short bursts of duration comparable with 
the life-time of electrons. On the other hand, independently 
on the assumptions, our results show that if particle accel- 
eration is mainly due to Alfven waves, then the existence 
of radio halos and HXR tails in massive clusters should be 
limited to achieve periods of < 10% of the Hubble time. 

vii) The temporal duration of the process of turbulent 
cascade, and the acceleration time scale of electrons are esti- 
mated to be about one order of magnitude shorter than the 
dynamical time-scale of a merger event. This implies that a 
temporal correlation is still expected between merging pro- 
cesses and the rise of the non-thermal phenomena in galaxy 
clusters. 
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